Near-thermal equilibrium with Tsallis distributions 
in heavy ion collisions 



J. Cleymans", G. Hamar 6 , P. Levai 6 , S. Wheaton a , 

a UCT-CERN Research Centre and Department of Physics, University of Cape 
Town, Rondebosch 7701, Cape Town, South Africa 

6 MTA KFKI RMKI, Research Insitute for Particle and Nuclear Physics, 
Po.B. 49, Budapest, 1525, Hungary 

£NJ , E-mail: jean . cleymans@uct . ac . za 

c3 



(N 



Abstract. Hadron yields in high energy heavy ion collisions have been fitted with 
thermal models using standard (extensive) statistical distributions. These models give 
| insight into the freeze-out conditions at varying beam energies and lead to a systematic 

consistent picture of freeze-out conditions at all beam energies. In this paper we 
' . investigate changes to this analysis when the statistical distributions are replaced by 

Cln' non-extensive Tsallis distributions for hadrons. We investigate the particle yields at 

■ SPS and RHIC energies and obtain better fits with smaller x 2 f° r the same hadron 

D | data, as applied earlier in the thermal fits for SPS energies but not for RHIC energies. 

<N 

i ■ 1. Introduction 

t> ' 

After many years of investigating hadron-hadron and heavy ion collisions, the study 
of hadron production remains an active and important field of research. The lack of 
detailed knowledge of the microscopic mechanisms has led to the use of many different 
models, often from completely opposite directions. Thermal models, based on statistical 
weights for produced hadrons [H El El H] , are very successful in describing particle yields 
at different beam energies [S1E1E1E], especially in heavy ion collisions. These models 
assume the formation of a system which is in thermal and chemical equilibrium in the 
hadronic phase and is characterised by a set of thermodynamic variables for the hadronic 
phase. The deconfined period of the time evolution dominated by quarks and gluons 
remains hidden: full equilibration generally washes out and destroys large amounts of 
information about the early deconfined phase. The success of statistical models implies 
the loss of such information, at least for certain properties, during hadronization. It is 
a basic question as to which ones survive the hadronization and behave as messengers 
from the early (quark dominated) stages, especially if these are strongly interacting 
stages. 

In the case of full thermal and chemical equilibrium, relativistic statistical 
distributions can be used, leading to exponential spectra for the transverse momentum 
distribution of hadrons. On the other hand, experimental data at SPS and RHIC 
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energies display non-exponential behaviours at high px- One explanation of this 
deviation is connected to the power-like hadron spectra obtained from perturbative 
QCD descriptions: the hadron yield from quark and gluon fragmentation overwhelms 
the thermal (exponential) hadron production. However, this overlap is not trivial. One 
can assume the appearance of near-thermal hadron distributions, which is similar to 
the thermal distribution at lower p?, but it has a non-exponential tail at higher pf. 
A stationary distribution of strongly interacting hadron gas in a finite volume can be 
characterized by such a distribution (or strongly interacting quark matter), which will 
hadronize into hadron matter. Tsallis distributions satisfy such criteria [TUl E] . In the 
next Section we will review the Tsallis distribution and emphasize the properties most 
relevant to particle yields. 



2. Tsallis Distribution for Particle Multiplicities. 

2.1. Relation between the Boltzmann and Tsallis distributions 

Neglecting quantum statistics, the entropy of a particle of species % is given by |9j 

Si = vJ ^(nf-nflnnf), (1) 

where the mean occupation numbers, nf , are given by 

nf(£)^exp(-^^). (2) 

with gi being the degeneracy factor of particle i. The total number of particles of species 
i is given by an integral over phase space of eq. (TSJ): 

Nt = v J Wr" m - (3) 

The transition to the Tsallis distribution makes use of the following substitutions [TU] 

ln(x) -> \n q (x) = — - , (4) 

exp(x) — > exp g (x) = [1 + (1 — g)^] 1 ^, (5) 
which leads to the standard result [TDl [EE] 



nf(E) = 9i 



l + (q-V El * 



9-1 

(6) 



T 

which is usually referred to as the Tsallis distribution [TQlfTT] . As these number densities 
are not normalized, we do not use the normalized g-probabilities which have been 
proposed in Ref. [UJ. In the limit where q — > 1 this becomes the Boltzmann distribution: 

\\mnJ(E)=n*(E). (7) 

q-*l 

The particle number is now given by 

Ni = v J ' Wf {n ' {E)) '- (8) 
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Figure 1. Comparison between the Boltzmann and Tsallis distributions. The figure 
on the left compares the two as a function of the energy E, keeping the temperature 
and chemical potential fixed, for various values of the Tsallis parameter q. The value 
q = 1.5 is the maximum value that still leads to a convergent integral in Eq. [8] 
chemical potential /i, keeping the temperature and the energy fixed. 



Note that q = 1.5 is the maximum value that still leads to a convergent integral in 
eq. (jHJ). A derivation of the Tsallis distribution, based on the Boltzmann equation, 
has been given in Ref. [12] . A comparison between the two distributions is shown in 
Fig. (00), where it can be seen that, at fixed values of T and /i, the Tsallis distribution 
is always larger than the Boltzmann one if q > 1. Taking into account the large pt 
results for particle production we will only consider this possibility in this paper. As 
a consequence, in order to keep the particle yields the same, the Tsallis distribution 
always leads to smaller values of T for the same set of particle yields. The dependence 
on the chemical potential is also illustrated on the right of Fig. [1] for a fixed temperature 
T and a fixed energy E. As one can see, the Tsallis distribution in this case increases 
with increasing q. The Tsallis distribution for quantum statistics has been considered 

in Ref. [HI [TH na ng. 

3. Relation between the Tsallis parameter q and temperature fluctuations 

The parameter q plays a central role in the Tsallis distribution and a physical 
interpretation is needed to appreciate its significance. To this end we follow the 
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analysis of Ref. [T7j and write the Tsallis distribution as a superposition of Boltzmann 
distributions 

- T(B) =r d (i) e - ( ^(i)- <9 » 

where the detailed form of the function / is given in [T7j. The parameter Tg is the 
standard temperature as it appears in the Boltzmann distribution. It is straightforward 
to show [T7] that the average value of 1/Tg is given by the Tsallis temperature: 

1 \ C°° . f 1\ f 1\ . f 1\ 1 



T B / ~ Jo d \Tb) \Tb) ! \Tb) = T' < 10 » 
while the fluctuation in the temperature is given by the deviation of the Tsallis parameter 
q from unity: 

1 I \ / 1 



Tb I I \T B 

q-l (11) 

A) 

which becomes zero in the Boltzmann limit. The above leads to the interpretation 
of the Tsallis distribution as a superposition of Boltzmann distributions with different 
temperatures. The average value of these (Boltzmann) temperatures is the temperature 
T appearing in the Tsallis distribution. This is the interpretation of the Tsallis 
temperature that we will follow. The other parameter in the Tsallis distribution, q, 
describes the spread around the average value of the (Boltzmann) temperature T. For 
q = 1 we have an exact Boltzmann distribution, for values of q which deviate from 1, 
we have a corresponding deviation. From this point of view the Tsallis distribution 
describes a distribution of (Boltzmann) temperatures. A deviation from q = 1 means 
that a spread of temperatures is needed instead of a single value. 



4. Thermal Fit Details 



In order to identify the energy dependence of the deviation from ideal gas behaviour, 
thermal fits were performed on yields measured at the CERN SPS in central Pb-Pb 
collisions at 40 AGeV, 80 AGeV and 158 AGeV (using the same data as analyzed in 
[H]) and yields measured at RHIC in central Au-Au collisions at \/s = 130 AGeV (using 
the same data as analyzed in [19]) and at y/s = 200 AGeV. 

In the CERN SPS fits, the thermal parameters T, /ig, 75 and R were fit to the 
data, while \lq and \is were fixed by the initial baryon-to-charge ratio and strangeness 
content in the colliding system, respectively. 

In the case of the RHIC analysis we again fit T, /jg, \x$ and 75 to the data. The 
use of mid-rapidity data here led to the relaxing of the constraints on fi$ and hq typical 
in analyzes of 47r data. Instead, fig was set to zero as justified by the observed 7r + /h~ 
ratio. The following expression was used to calculate primordial particle yields, 

nf(£,7 s )=7f 



E 



T 



-j-m— -v 



(12) 
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Figure 2. The x 2 /d.o.f. of the fits as a function of the Tsallis parameter q. 



where \Si \ is the number of valence strange quarks and anti-quarks in species i. The value 
7 S = 1 obviously corresponds to complete strangeness equilibration. All calculations 
were done using the THERMUS package [18J. 

5. Results and Conclusions 

The most surprising result of our analysis is shown in Fig. (T5]): the quality of the fits, 
as measured by the x 2 /d.o.f., improves at first as the Tsallis parameter q increases. It 
reaches a minimum value around q ~ 1.07 for SPS beam energy of 158 AGeV. This 
behaviour is repeated at other SPS energies with the minima at slightly different values 
of q, i.e. 1.08 for 80 AGeV and about 1.05 at 40 AGeV beam energy. This behaviour 
is not seen at RHIC energies. Clearly, changes in the Tsallis parameter q have only a 
small negligible effect on the y 2 values at RHIC energies, of course, this still leaves open 
the possibility for q values larger than 1 [20J. However on the SPS data the effect is 
substantial and changes the interpretation substantially. One possible interpretation is 
that at SPS energies fluctuations in the freeze-out temperature are substantial. 

Recently [21] a coalescence model with a Tsallis distribution for quarks was used 
to fit the transverse momenta spectra measured at RHIC. This fit does not include 
decays from resonances and therefore cannot be compared directly to ours since decays 
can substantially modify the transverse momenta, also the emerging hadrons are not 
in a Tsallis-type equilibrium gas, which is an assumption of the present analysis. The 
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Figure 3. The chemical freeze-out temperature (left figure) and the baryon chemical 
potential (right figure) as a function of the Tsallis parameter q. 



authors obtained values for the Tsallis parameter q which are remarkably similar for all 
particle species considered, i.e. q ~ 1.2 which cannot be excluded by our analysis. 

The freeze-out temperature T decreases, as expected, with increasing values of q. 
This can be understood from the fact that the Tsallis distribution is always larger than 
the Boltzmann one (as long as q > 1). Hence, in order to match the same particle yields 
one has to adjust T to lower values. This is seen at all energies in Fig. [3j However, the 
drop in T, turns out to be quite drastic numerically. In fact, the decrease in particle 
numbers has to be compensated by increases in all other thermodynamic variables. The 
(modest) increase in the baryon chemical potential is shown in Fig. [3] on the right hand 
side. 

The strangeness non-equilibrium factor 7 S as shown in Fig. HI It is interesting to 
note that the Tsallis distribution leads to a much better chemical equilibrium than the 
corresponding Boltzmann distribution with q = 1. In all cases considered the 7 S is very 
close to the chemical equilibrium value of 1. Clearly, the use of the Tsallis distribution 
in relativistic heavy ion collisions calls for a reevaluation of the understanding gained 
from previous analyses El E] • 



Near-thermal equilibrium with Tsallis distributions in heavy ion collisions 
1.3 



1.2 



1.1 



0.6 



0.5 



1 I 1 I 1 

SPS 158 A GeV 

SPS 80 A GeV 

SPS 40 A GeV 

RHIC 130 A GeV 



T 



RHIC 200 A GeV 




_L 



J_ 



1.02 1.04 1.06 1.08 
Tsallis Parameter q 



1.1 



1.12 



Figure 4. The strangeness non-equilibrium factor 7 S as a function of the Tsallis 
parameter q. 
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